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In this reply we answer a comment by A. Dhar on our 
Letter |l[ . In that paper we studied thermal conductivity 
in a one-dimensional gas of hard point alternating-masses 
particles, demonstrating that Fourier's law holds both in 
its statical and dynamical aspects. 

Dhar's main point is about the asymptotic behavior 
of the total energy current self correlation function C{t). 
He fits C(t) ~ t~°- 83 , implying an infinite thermal con- 
ductivity while we obtained C(t) ~ t~ 13 and thus a 
finite k. In our opinion, Dhar's result is a consequence 
of autocorrelations due to finite size effects. In order to 
show this point, we plot in fig. [jja C(t) for a system size 
N = 1000. Here we can study two different regions: (1) 
one for ln(i/*o) G [8,9], and (2) other with ]n(t/t ) > 10. 
A power law fit to the first region yields C(t) ~ 
while a fit to the second region yields C(t) ~ i~ 88 , 
very similar to Dhar's result. In fact, similar behav- 
ior has been measured by Savin et al0. However, we 
think that only region (1) corresponds to the infinite sys- 
tem asymptotic behavior. As an example of the pre- 
vious statement, we can study the asymptotic behavior 
of the local energy current self correlation function c(t) 
for the equal masses gas. Following Jepsen it can 
be shown analitically, after a lengthly calculation, that 



(1 ) : ln|0(t)|=1 2.6(0.1 )-2.98(0.02)ln(t/t ) 

(2) : ln|c(t)|=-0.1 (0.1)-0.83(0.02)ln(t/t ) 


ln(t/t ) 


FIG. 1: (a) C(t) for the different masses system with N — 
1000. (6) c(t) for the equal masses system with N = 500. The 
insets show the results of a power law fit for both regions (1) 
and (2). to is the mean collision time. 


c(i) ~ t~ 3 for this system. We have measured c(i) jjj for 
equal masses in a finite system in the canonical ensemble 
(not in the zero-momentum ensemble, contrary to Dhar's 
comment). Fig. |[b shows c(t) for N = 500. It is remark- 
able that we can also define here two different regions: (1) 
one for ln(t/to) € [5.1,5.8], where a power law fit yields 
\c(t)\ ~ t~ 3 , and (2) one for ln(t/t ) > 6, where a power 
law fit yields \c(t)\ ~ t~ 0S3 . We recover the theoretically 
predicted asymptotic bulk behavior in region (1), while 
region (2) should be due to finite size effects. Moreover, 
it is intringuing that the finite size time decay exponent 
(~ 0.83) is almost the same both in the different masses 
case and the equal masses one. This fact points out the 
existence of an underlying common finite size mechanism, 
responsible of this spureous long time decay In conclu- 
sion, coming back to the different masses case, we think 
that the above example indicates that only region (1) of 
fig. [l].a represents the asymptotic bulk behavior. Hence, 
any conclusion about system's conductivity derived from 
region (2) should be misleading. 

Let's clarify now some other minor points raised in 
Dhar's comment. First, we do not observe linear tem- 
perature profiles. They are linear in the central region 
and curved near the boundaries, which is consistent with 
the finding of finite conductivity. Q Second, the validity 
of our deterministic heat bath and our local tempera- 
ture measure has been carefully tested. On the other 
hand, it can be shown that C(t) ~ Nc(t) + J2i^i c m(0- 
Hence, for regular systems, where non-local time corre- 
lation functions Cij(t) decay fast enough with distance, 
one expects a similar long time decay for both C(t) and 
c(i). However, there are anomalous systems, as the (non- 
ergodic) equal masses gas, for which c,,;(t) decays very 
slowly, or does not decay at all, and thus C(t) and c(t) 
behave completely different. 

In conclusion, we firmly confirm, after a global, con- 
sistent analysis of the problem, our previous results pj, 
i.e. that our one-dimensional system has a finite thermal 
conductivity in the Thermodynamic Limit, thus obeying 
Fourier's law. 
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